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INTRODUCTION 

Lithospheric plates descend into the mantle at deep ocean trenches. TItis observation is a 
fundamental part of the theory of plate tectonics. During their descent these slabs undergo an 
evolution in thermodynamic state which, presumably, includes density changes. It is of prime 
geophysical importance to define the anomalous density structure of the slabs. 

Sea-going gravimetric measurements have been used by quite a number of workers (e.g„ 
Segawa and Tomoda, 1976; Watts and Taiwan!, 1975; see also Grow & Bowin, 1975) to study 
density structure, including that of slabs, in the vicinity of subduction zones. These studies have 
had limited success in resolving slab density structure. Gravity anomaly profiles across trench/island 
arc complexes are characterized by a high and an adjacent low, each of large magnitude (typically 
~200 mgal) and short wavelength (~200 to 300 km). As Watts and Taiwan! (1975) indicated, 
these gravity anomaly signals are so dominated by crustal effects that it is difficult or impossible 
to recognize the gravitational contribution of the slab. One expects the slabs to have a contribution 
of generally longer wavelength, by reason of tlieir depth, than that of crustal effects. Fortunately 
direct observations of sea surface topography (effectively the geoid for these purposes) which are 
provided by Geos-3 and/or Seasat altimeters afford a more precise resolution of the gravity field’s 
intermediate wavelength (600-4000 km) component (Cliapman and Taiwan!, 1979). Furthermore, 
shorter wavelengths due primarily to near surface (e.g., crustal) effects arc attenuated in the geoid 
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representation. The gravitational effects of slabs should be primarily of intermediate wavelengths 
and longer. Therefore the altimeter data are particularly useful in defining the density structure of 
slabs and in elucidating the global correlation between subduction zones and gravity highs. This cor- 
relation has been pointed out by Kaula (1972), McKenzie (1969) and Chase (1979) among others. 
Note that in this study sea surface ; npography will be assumed to represent geoid heiglits. This is 
due to the fact that oceanographic effects, constituting the difference between the geoid and the sea 
surface, generally should not exceed 1 or 2 meters. 

A major difficulty, and in certain cases an asset, with interpreting geoid anomalies arises due to 
the fact that the geoid height at any one point is influenced by density variations over a relatively 
extensive region. Gravity anomalies on the other hand are largely due to nearby crustal density 
effects as indicated by the general utility of the simple Bouguer reduction. Therefore in interpreting 
the geoid in a given region, such as the southwest Pacific, it may be important to include slabs of the 
adjacent region (e.g., the New Guinea— New Britain— Solomon slabs, see below). The geoid in a given 
region may be particularly responsive to somewhat distant density anomalies which are not locally 
compensated. 

On the positive side, longer wavelength density anomalies more clearly manifest themselves in 
the geoid. This may be seen by examining Stokes’ integral formula in surface spherical harmonic 
form (Heiskanen and Moritz, 1967). Specifically, subducting slabs (long along strike) become sub- 
stantial, perhaps dominant, contributors to geoid undulations as a direct result of their length. If 
one assumes that slabs 'Vi orldwide possess a net linear mass anomaly (i.e,, are not regionally compen- 
sated) of the order 10*^ gm/cm, it can be seen in what follows below that slabs contribute, down to 
quite low degree harmonics, to the global geopotential power spectrum. In other words, slabs may 
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be an important source of large long wavelength gravity anomalies. As a counterexample, the South 
Sandwich slab is so short-along strike, that its effect on the geoid is not very noticeable. 

The followhl,® are some major simplifying assumptions that have been invoked: 

• The Benioff zone (zone of seismicity) completely defines the position of the anomalously 
dense litho.spheric slab. 

• Lithospheric slab thickness is not significant for these purposes and an individual slab may 
be considered to be condensed onto an infinitesimally thin, median surface of uniform 
mass per unit area. 

• The contribution of crustal density variation may to a reasonable approximation be ignored. 

• Thermal litliosplverc associated with young marginal basins in the southwest Pacific and 
with the Pacific plate proper will be considered to be isostatically compensated and to have 
tlie same geoidal contribution as comparably aged normal oceanic lithosphere elsewhere 
(see Haxby and Turcotte, 1978, also see below). 

A number of workers (including Toksoz et al,, 1971, Griggs, 1972, Mincar and Toksoz, 1972, 
Scliubert et al., 1975) have developed models of the thermal behaviour of the lithosphere as it sub- 
ducts into the mantle, Tlieir results indicate that the effects of thermal contraction and elevation of 
mineralogical phase transitions should produce positive density anomalies in tlie slab. Tliese density 
anomalies should vary substantially with depth, In the models that follow, density anomalies are 
assumed uniform with depth; discrepancies between model results and observations do not readily 
indicate just how density anomalies do vary with depth. Nor has it been possible to resolve any 
regional compensation of tlic slabs such as Griggs (1972) proposed. 
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In the case of the southwest Pacific, it is important to consider that different slabs liavc dif- 
ferent overall density anomalies. This has been done. The model suggests that the New Hebrides 
slab possesses an overall average density anomaly of +3,2 X 10® gm cm*^, which is about three 
times the value obtained for the Tonga-Kermadcc slab. Of course a unique inversion or density 
determination is not theoretically possible. 


LINE SOURCE MODEL-GREAT CIRCLE SEGMENTS 


As a prelude to the detailed models which follow, let us consider a rather idealized model of 
the global gravitational effect due to a slab of substantial length along strike, Consider the slab to 
be a line source of density, Pj, mass per unit length. Also consider it to lie along a great circle seg- 
ment of arc length 2ri//g at some depth, d, (cf., Figure 1) in a spherical Earth. The change in gcoid 
height at any point r' is given by 


6N (t') « 



di// 


—♦-4 

r - r 


( 1 ) 


where r = | r |, G is the universal gravitational constant and g is gravitational acceleration at sea level. 
The maximum geoid height will occur over the center of the line source at r' = r'x (i.e., >//', = 01. 

For this case equation (1) may be rewritten as, 


5N(r',^' = 0) 


where 



( 2 ) 


and the substitution <f> = ( 7 T - i //)/2 is used. The above integral is simply the difference between a com- 
plete elliptic integral of the first kind and an incomplete one. A general evaluation of potential 



Figure I . Geometry of a line source having the shape of a great 
circle segment. 


integrals of elliptic form can be found in Morrison (1980). For shallow depths (d <^r') equation (2) 
may be surprisingly well approximated by the result for a straight line source (length 2r\i/Q): 

2Gpj / 

6N (r', 0) = In IVl + + f 

g \ 

_ '•'‘'o 

where ? = "T“ 



This result is easier to evaluate. For a line source at 300 km deptli, tlie maximum geoid height, 8N, 
computed using the straight line approximation (2) deviates by less titan 3 percent from the actual 
value of equation (1). These two results as functions of source length—are shown in normalized 
form in Figure 2a. 
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The corresponding gravity anomalies (actually distrubanccs), due to this line source are given 
(to an excellent approximation) by tlie straiglit line result: 

2p^G 

5g.x(f^0) (4) 

dVl +(dlt\pQ)^ 

where d « r. Tills expression is evaluated in nonnalized form in Figure 2b for the depth 
d = 300 km which was used in computing Figure 2a. 

It is clear from Figure 2b why the two-dimensional approximation (i.e., that the structure is 
of infinite extent in the direction normal to the cross-section) is f,o use ul in modelling gravity 
anomaly profiles which transect linear features. However, this same two-dimensional approxima- 
tion (cf, Figure 2a) is not appropriate for modellint geoid undulations. The length of the structure 
(in the direction normal to the cross section) has a substantial influence on the magnituds of the 
geoid anomaly which it produces. For example, the line source model suggests that the Tonga- 
Kermadec slab (rouglily corresponding to sliould liave a geoid anomaly amplitude about 

three times that of the South Sandwich slab (V'q'^^®) assuming the same slab depth (actually the 
South Sandwich appears to be shallower) and linear mass density, Pj . The nominal wavelengtli of 
this larger geoid anomaly will also be correspondingly larger. In contrast to the computed geoid 
anomalies, tliese two line source models should produce gravity anomaly profiles of about the same 
magnitude and shape. This is a simple realization of what was pointed out in tlie Introduction; 
longer wavelength mass anomalies manifest tliemselves more strongly than short wavelength anom- 
alies in the geoid representation. 
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Let us now consider tlic case of a complete great circle source (ring source); i.e., ** 1 80” in 

Figure 1 . The axis or pole of this ring will be located at an arbitrary longitude, X', and colatitude 
O', If O ' « 0, for example, the ring will be equatorial. Using a derivation similar to that used by 
Pollack (1973) for point and cap masses, the potential due to tlic ring source can be represented as 
a spherical hannonic expansion. This potential is given at an arbitrary observation point (r, X, 0) 
by: 


AVJr.X.O)*- 


G 2-n PgUp P| 


n"0 


(o)/RV‘r 1 

7T\ \T) + ^m«^'.^')Smn(0.^)J(5a) 


where R^, is the radius of the ring, and are the fully nonnalized surface spherical har- 
monic functions of degree n and order m, and P„ (O) are the unnormalized Legendre polynomials 


evaluated at cos 7r/2, Excepting G, pg, and R^, , all symbols in (5a) conform to the definition used in 
Heiskanen and Moritz (1967). 


Let us assign to this ring source a linear density, pg, which is comparable to expected density, 
i.c., mass per length along strike of a typical slab. Then a power spectrum can be computed and 
compared to an empirical measure (e.g., Kaula’s rule of thumb) of the earth’s geopotential power. 
The earth’s anomalous gravitational potential, AV, may be expressed as: 


AV 




■^iim ^nm 


+ »^n.n S„„, (0 


.X)j 


(5b) 


where M and R are respectively the mass and the nominal radius of the earth. Alsol,,,,, and 
are the familiar coefficients of this normalized expansion. The degree variances for this geo- 


potential are given by: 



(5c) 
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An empirical law due to Kaula (1968) (cf, Lambeck, 1976) provides a rough estimate of for 
harmonics up to degree 20, according to tlie law: 






o/ ~ (2n + l)10*"’/n’ 


(5d) 


Figure 3 includes a plot of root variance, given by this law. Comparable degree variances for 
the ring source can be determined. Excluding degree ri * 0, 1, the expression (5a) for the potential 
AVq can be written exactly in the form of (5b) by defining the coefficients ® as follows; 

/^nm \, . ^ffPgRp /R^" /^im > 


Then using (5c) ring source degree variances can be computed and directly compared to 
tlic observed variances given approximately by (5d), This comparison is sliown in Figure 3, Tlie 
root degree variances for ring sources with two differing linear densities, pj, are shown. Note in 
computing degree variances from (5e), one can either set 6' = 0 or take advantage of the orthog- 
onality property of both (0', X') and S,„„ (0', X'). The ring source is placed at depth 
d = 300 km; i.e,, s 6070 km where a nominal earth radius R - 6370 km is assumed. Also note 
that odd-degree variances of the ring source potential are zero. 


A noteworthy aspect of Figure 3 is the near congruency of one ring source variance curve 
(Pj = 2 X I0‘^ gm cm’* ) to that of observed or empirical variance curve, interestingly, the mass 
per length, which was determined below, to best fit the New Hebrides slab, works out to be approxi- 
mately 2 X 10*^ gm cm’* . It is not clear that this density is typical of slabs worldwide, nor is it 
clear that this ring source model can provide a rougli estimate of the power spectrum due to the 
global slab distribution. Overlooking these problems one miglit interpret this near congruency of 
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DEGREE ,n 

Figure 3, Plots of root degree variance as a function of degree for; i)Kaula’s rule-of> 
thumb equivalent; (ii) a ring source of density - 1 X 10*^ gm cm'* and 

depth 300 km; (iii) a ring source of density Pj = 2 X 10*® gm cm'* and depth 
300 km. 
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spectrum curves as an indication that slabs contribute substantial power to the observed geopotential 

degrees ^ 20 ), 


SURFACE SOURCE MODEL 

The slab is assumed to be an inrmitesimally thin surface. This assumption ic reasonably jus- 
tified, Upward continuation of gravitational potential from slab depth to sea level washes out 
virtually all the manifestations of tlie slab's finite thickncss-except perhaps those due to the slab’s 
shallowest portions. 


A surface, S, of unspecified analytical fonn, can be fit to the seismic data defining the Beni- 
off zone. This slab surface will produce a geoid anomaly at T' given by: 




// 


gds 

fT'-T" 


(6) 


where ds is an element of surface area, It is the position vector of a surface element, and a is the 
mass per unit area. Let us specify that this surface, S, be polyhedral and be composed of m planar 
triangular facets, S,, each having a uniform mass per unit area a. Then equation (6) becomes: 


«N(x',y',z') « 


or 


■ ^ j ‘ (y' - y)* (*' - ds 

i»l * ( s, 

* m 

5N(x’.y'.z’) - 6N, 




(7a) 


(7b) 


By defining the apices of all facets (see Figures 4 and 5), one effectively defines the entire slab 
surface. 
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Figure 4. Geoinetiy of a typical planar triangular facet showing transformed (u, v, wj arid untransformed 

(x, y, 2) coordinates systems. 


The integral expression' for eacli facet’s contribution, 5Nj, can be evaluated using methods 
applied by Chapman (1979), Barnett (1976) and Paul (1974), The integration in equation (7) can 
be accomplished analytically after a coordinate transformation from an arbitrary set (x, y, z) to 
one (u, V, w) in whicli the facet is coplanar with u, v (see Figure 4). This transformation is described 
by Paul (1974) and Chapman (1979); it consists of a rotation and a translation and is given by: 




^"11 "12 



V 

■ 1 

1 

iij, iijj nj3 


y-vo 



i”31 '^32 *’ 33 ) 




wiiere (x^ , y,_, , ) is the position of the observation point and n|j is the familiar direction cosine 

matrix, For example, n. 2 is y component of the unit vector ft, . Note from Figure 4 that ftj = 
ft, X ^2 ftj = 1I3 X ft, . Tenn 5N, of (7b) becomes-in transformed coordinates- 


5N. 


/ 


du dv 


AS, 


[u* +v® + w®)''^' 


(9) 
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Figure 5, Contours defining the position of the slabs, Siabs are considered to be infinitesimally thin. Crosses 
define apices of individual blocks, each of which comprises two triangular facets. Finely dashed 
line defines slabs (segments) which correspond to precisely located Benioff zones. See text for 
etails. 
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By applying Stokes’ theorem to the Integral In equation (9), it can be rewritten as: 

Go. 


6N. 


-5- J " fin |*u + + w’ J 


dv 


( 10 ) 


where C is the perimeter closed contour and integration is performed In a riglit-handed positive w 
sense. For the triangular facet sliown in Figure 4, one can integrate (10) and obtain: 


Go, 


5N. 


( h fin r«,P,(a^l)-‘''-K«’ + l)''V 3 +R 3 “j 

) [uj + R, J %l + it Pj ia\ + 1)-^ + ia] + l)’-^ V, + RjJ 




(a\ + 1)''^ 


fin 


Us h H + 


iy ’''‘+(«5 + i )'''' v , + R ,' 


*3 '•“3 i) ^ (<*3 '^3 B3 


+ W 


+ tan' 


(«, - /3, V 


tan' 


2'3\ . 


I- tan'' 


w R, 


(• 


wR, 


\ W R. j 


tan' 


f 


wR, 


( 11 ) 


where (u, , v, , w), (u^, Vj , w), and (Ug, Vj, w) are the locations of apex 1, 2, and 3, respectively 
(note Vg = V, and Wg = w^ = Wj = w) and where: 

«j = (Ug - Ug/(Vg - V,) , 

“3 = CUj - Ug)/(V, - Vg), 

p. — U| + ttj Vj , i — 2, 3, 

R. = (Uj^ + Vj* + i = 1,2,3, 


This result (1 1) is similar to that obtained by Chapman (1979) and Barnett (1976) for a three- 
dimensional body; i:he only difference being that this result is developed for a polyhedral surface 
source. 
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Sutftices representing the Tonga-Kennatlcc slab and the New Hebrides slab are shown in 
Figure 5, The Tonga«Kcrinadec slab is blocked in 0.5° latitude intervals along each pair of adjacent 
depth contours. The contours are adapted from those of Billington (1979). Each block is com* 
posed of two triangular facets. The facets on a pair of typical blocks are shown in Figure 5, All 
other blocks are sinularly facetted; this is not shown. There are a total of 944 triangular facets in 
this slab. Tlie shallow portion (<310 km depth) of the New Hebrides slab is blocked in 1.0° 
latitude intervals and totals 68 triangular facets. These depth contours arc adapted from those of 
Isacks and Barazangi (1977), 

The depth contours of Billington (1979) and Isacks and Barazangi (1977) define the position 
of the top of the Benioff zone. In our slab model (sec Figure 5), these seismicity contours have 
been shifted using a slab half-thickness of 25 km and a single overall average slab dip which is 
applied to all contour excepting the trench axis contour. Tlie trench axis contour is shifted directly 
downward. 

It is important to include continuations of the Benioff zone westward from the northern 
New Hebrides and southward tiom the southern Kermadec. These zones (slabs) are relatively dis- 
tant from the points at which their geoid contribution is to be computed. Over such distances, 
short wavelength gravitational effects are washed out and only crude models (coarse blocking) arc 
needed. The Benioff zone(s) of the Solomon-New Britain-New Guinea region has been crudely 
represented using cither one or two blocks (i.e., one or two depth intcivals) for each 2° longitude 
interval from I65°E westward to 135°E (see Figure 5), These two-dimensional blocks represent 
approximate fits to the seismic zone determinations of Pascal (1979) and Johnson and Molnar 
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(1972). Sunilarly, a crude fk to the seismicity distributions derived by Hamilton and Gale (1968) 
for North Island, New Zealand, provide an extension of the southern Kennadec slab to 41.5®$ 
latitude. A portion of this fit is shown in Figure 5 (hlplwrequency-dasli contours). 

A rectcji)*ular contour at 650 km depth under tlie northern Fiji plateau encloses three blocks 
which crudely fit the seismic zone detailed by the Isacks and Barazangi (1977) and Pascal ct al 
(1978). Both of these papers suggest that this zone may coincide with a remnant segment of sub- 
ducted Pacific plate. The Vitiaz trench is thought to be a legacy of this extinct subduction zone. 
This rectangular segment will be lumped with the Tonga-Kermadec slab when assigning density 
anomalies and computing geoids. 

The 550 km contours under the soutliem Fiji Plateau define a trapezoidal slab segment which 
roughly fits a seismic zone detailed by Barazangi and Isacks. Both they and Pascal et al (1978) 
suggest that this zone probably is related to ongoing New Hebrides subduction. Therefore this 
550 km deep segment will be lumped with the shallower New Hebrides slab in assigning densities 
and computing geoids. Evidence that there is indeed a discontinuity (as idealized in Figure 5) 
between the shallow and deep slab portions is presented by Barazangi st al (1973), 

Theoretical geoid anomalies computed from this slab model are shown in Figures 9-11. Com- 
parisons with observations are contained in the Results section. It can be seen tliat the character 
in the computed geoid anomalies is largely attributable to the Tonga-Kermadec and shallow New 
Hebrides slabs. In other words, the more precisely located slabs (Benioff zones), not the crudely 
located slabs, provide the computed geoid with its important shape and character. 


A MODEL OF THERMAL LITHOSPHERE ASSOCIATED WITH MARGINAL BASINS 

Recent studies (Kaxby and Turcottc, 1978;Haxby, 1980; Crougli, 1979) have shown that 
observed gcoid lieiglits tend to vary systematically with the age of normal oceanic lithosphere. 
In particular, Haxby and Turcotte (1978) sliowed that observed geoid lieights tend to decrease 
linearly with increasing lithospheric age out to 80 m.y.B.P. in good agreement with the half-space 
simple conductive cooling model, Crough (1979) showed that there are significant gcoid height 
changes across fracture zones which separate lithospheric segments of differing ages, A similar 
age contrast exists across deep ocean trenches, It appears that for a number of trenches, geoid 
heights increase as the observer crosses from areas seaward of the trench to those landward of the 
trench. Part of this height increase is due to the contrast between the young, thin lithosphere of 
the marginal basins ami the typically older, thicker oceanic lithosphere seaward of the trench. 


It will be assumed a priori, that thermal lithosphere associated with four marginal basins 
(sec Figure 6) of the soirthwest Pacific can be modelled as normal oceanic lithosphere. More 
exactly, the following relations apply: the depth-age relation from thermal isostasy: 


w « \Vn + 


2p a(T -T)/xt\^ 
^ m o' I 1 

(p„rPw‘) 


( 12 ) 


due to Davis and Lister (1974), among others, and the associated geoid height-age relation 


27rG 

AN = - p 

• II 


2p a(T -T )’ 

^ K O' 


(13) 


due to Haxby and Turcotte (1978). In the two preceding equations, T,,^ is the initial isotliermal 
mantle temperature; T^ is the invariant surface temperature; is mantle density at temperature 
seawater density; k is the mantle thermal diffusivity; a is the volumetric coefficient of 
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Figure 6. Dashed perimeters define the extent of four marginal basins. These basin definitions are used in 

constructing the two-layer doublet shell model of the thermal lithosphere. See text and Table 1 for 









thermal expansion; w,^ is the depth to the right crest; t is the particular lithospheric age; and t,^ is 
age of lithosphere used as a reference (i.e„ corresponds to zero geoid height). 

Marginal basins differ distinctly in character from normal oceanic lithosphere; althougli they 
both evolve via crustal extension. Certainly the assumption of similarity, outlined above, is some- 
what tenuous. For an individual marginal basin, it is quite difficult to observe the depth-root age 
relation (12), This difficulty is perhaps due to characteristically rough topography (cf„ Karig 
et al., 1978), However, the South Fiji Basin, the Lau-Havre Basin, and the Fiji Plateau (see 
Figure 6) each have paired values of average depth-average age which fall nearly on the depth-age 
curve for the North Atlantic or North Pacific (see depth anomalies in Table 1), Other marginal 
basins, such as the Mariana Trough, lying outside this study’s region of interest, are substantially 
deeper, on average, than comparably-aged normal oceanic lithosphere (Karig et al., 1978). The 
rough agreement between the normal oceanic lithosphere’s age depth relation and that of the 
marginal basins in this study provides justification for applying the model given in (13), 

In computing a map of geoid anomalies due to marginal basin/oceanic lithosphere contrast, 
all geoid heights refer to a zero datum which corresponds to the Pacific plate seaward of the Tonga- 
Kermadec trench. In other words, the Pacific plate reference age, t^., in (13) is 75 m.y.B.P, 

Equation (13) will not be used explicitly to compute geoid anomalies. This expression is an approx- 
imation which is valid for long wavelength features (> 200 km). 


As Haxby and Turcotte (1978) have indicated, it can be expressed as 


2rrG 

AN = M 

g 


(14) 
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Table 1. 

Marginal Basins of the Southwest Pacific 


Busin 

♦Period 
of Activity 
(10® yrs) 

♦♦Mean 

Depth-W 

(km) 

Representative 
Age-t 
(10® yrs) 

Dipole 
Moment, M 
(gm cm’‘ ) 

Depth , 

Anomaly 

(km) 

Fiji Plateau 

0-3‘/0-9Bp2 

2.9® 

2,016,0 

2.73 X 10^^ 

0.0/+0.3 

South Fiji 

28-35 BP^ 

4.6’ 

31.5 

1.63 X 10*2 

-0.1 

Lau Havre 

0-3,5 BP^ 

2,6® 

2,3 

2.71 X 10*2 

+0.3 

New Hebrides 

45-52 BP® 

<4.9 

48.5 

0.99 X 10*2 

- 


♦Periods of activity from'. 1. Chcrkis et al (1978); 2. Falvey (1975); 3. Watts and Weissel (1977); 4. Weissel (1977); 
5.Welsselctal(1977). 

♦♦Mean depths from: 6. and 7.Karig(1971);8. Toksozand Bird (1977), 


where fj. is the first moment, in depth, of the density distribution. Comparison with (13) yields the 
expression for moment 


= -T„) 


. 2p,„«(T„.Tn 

« 1 f (t^. 

L J 


t) 


(15) 


Tlie geoid anomaly map (e.g,. Figure 7) is computed by constructing an equivalent doublet shell of 
moment, //, as given above. The moment for each of the four basins is listed in Table 1 , The dipole 
moment can be decomposed into a source* term, a, given by 


» = 20n,“(T„, (1“-.“) 

and a separation distance d given by 

r 2p.„a(T,-T^)1 (t,-t) ) 

2 L (P^-Pjn J(t/‘-t'/') 


20 





sucli tlitH n * od, Two polyhedral surfaces are used to fonn the doublet shell. For each of the four 
basins in Table 1, the upper surface of areal density +o (see equation 16a) is placed at mean depth 
w", The lower surface of density *o is placed at depth w + d, where d is given by I6b. The gcoid 
anomaly map (Figure 7) is then computed using the facet integration (1® X 1® blocking) and sum- 
mation described in the preceding section. 

This technique results in a smooth geoid even at the edges of the basins. If equation (13) 
were applied directly, the resultant geoid would have discontinuities at the edge of the basin, 

The geoid computed from the doublet shell does not contain the spurious high frequency compo- 
nents which will occur in a geoid computed directly from (13). 

In calculating the dipole moment, separation distance, d, theoretical deptli, w, etc, for each 
basin a representative age (t in Table 1) is used. Also Parsons and Sclater’s (1977) best fitting 
parameters (North Pacific determinations) were used. They arei -3.33gmcnr®, fc = 8.0 X 
10'^ cm*,s'‘ , T„, = 1333° and 3.3 X 10‘® °C . The representative age is derived from 

published studies of the particular basin’s magnetically expressed history (see Table 1 ). It is not 
a mean age. It is the age associated by the age-depth relation (12) with the predicted mean depth. 
The predicted mean depth is computed using (12) and assuming a uniform spreading rate. Note 
that actual mean depths are shown in Table 1, These depths are due to Karig (1971) and Toksoz 
and Bird (1977). A mean depth for the New Hebrides basin was not included, as this basin’s 
apparent complexity made a meaningful determination difficult. 

Note that map coverage in Figure 7 is not complete. There are a number of different basins 
in this 40° X 40° region which have not been sufficiently studied to obtain estimates of the thermal 
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lithosphere’s age. Therefore# these regions (e.g. the southwestern corner) appear as gaps-l.e., 
they are shown as areas of zero geoid anomaly or lithospheric age contrast, 

Note also that the effect of a young basin on geoid heiglits is clearly localized. This results 
from the assumption of an isostatically compensated thermal lithosphere. 

RESULTS 

Figure 8 contains a map of sea surface topography In the southwest Pacific derived from 
Seasat-I altimeter data. This map was constructed from data obtained between July 28 and August 
1 5, 1978 and is due to J. G. Marsh et al (1980). Heiglvts are referenced to an ellipsoid of flattening, 
F - 1/298.255, and semi-major axis, a « 0,378,140 meters. Small crossover lieight differences 
permitted direct contouring of the along-track data into map form. Tides were not explicitly 
removed in obtaining the heights in this map. The contribution of oceanographic effects to this 
representation of sea surface topography can be expected to be generally less than one meter in 
magnitude. For the purposes of tliis study, the map (Figure 8) may be considered to represent 
the marine geoid. Wavelengths sliorter tlian 200 km are effectively filtered out of this map, 

Figure 9 represents the computed geoid anomaly due to the surface source model (as 
previously described) of the regional slab distribution. The assumed density anomalies are: (i) for 
the Tonga-Kerinadec slab, 1,0 X 10® gm cm"* and (il) for the New Hebrides slab, 3.2 X 10® gm 
cm** . fn this case, the New Guinea-New Britain-Solomon (NG-NB-S) slabs of the adjacent region 
are assigned the same density as theTonga-Kermadec, As explained earlier, deep slab segments 
under the northern and southern Fiji Plateau are lumped, for density assignment purposes, witli 
the Tonga and New Hebrides slabs respectively. 
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Figure 9. Computed geold hciglit contours in meters due to the slab model. Densities arc: 

1.0 X 10® gm cm'® for the Tonga-Kcrmadcc and NG-NB-S slabs; 3.2 X 10® gm cm*® for 
tlic New Hebrides slab. Dashed lines arc ground tracks of satellite passes used in Figure 
12. Dotted lines sliow the positions of the Tonga-Kcrmadcc and New Hebrides trench 



Ill comparing these regional computed geold anomaUes (Figure 9) with the observed sea sur- 
face (Figure 8), agreement In shape and character can be easily seen. Note in particular that the 
twin local maxima at longitude '^169®E, latitude -^IS® and longitude -184®E, latitude ~22®. 

As will be seen in Figures I2a*12e, one obtains a much better agreement by including the geoidal 
effect of litliospheric age contrast between young marginal basin areas and those of older Pacific 
lithosphere. Furthermore, this model ignores crustal and bathymetric effects. To some degree, 
such effects are represented in the sea surface map, For example, weak lows tend to be positioned 
over the Tonga-Kcnnadec and New Hebrides trench axis. 

The significant character in the computed geold (Figure 9) is largely attributable to those 
slabs (slab segments) which possess more precise fits, These fits result from Benioff zone determi- 
nations of Billington (1979) and Isacks and Barazangi (1977);they arc described in Figure 5. The 
geoidal effect of these slabs is shown in Figure 10, Clearly the twin maxima in the preceding com- 
puted gcoid (Figure 9) aiisc due to these slabs. 

Figure 1 1 shows another computed regional gcoid. This computation differs from that of 
Figure 9 only in that the density of the New Guinea-New Britain-Solomon (NG-NB-S) slab complex 
has been raised from 1.0 to 2,6 X 10* gm cm‘^ and the density of the New Hebrides slab lowered 
from 3,2 to 2,6 X 10* gm cm'® . Tliis gcoid is comparable to that of Figure 9 except that there is 
less closure of the contour lines in the northwest comer of this map (Figure 1 1). In this particular 
aspect the gcoid of Figure 1 1 provides sliglitly better agreement with the observed geoid. In fact, 
global sea surface maps or gcoids display a rather pronounced high over the New Guinea area. This 
suggests that the choice of a higlier (2.6 X 10* gm cm’®) nominal density anomaly for the NG- 
NB-S slabs may be closer to correct. However, the geoids computed in Figures 9 and 1 i arc not 
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sensitive to tlie particular density anomaly selected for these slabs. To obtain a reasonable density 
anomaly estimate for the NG*NB-S slabs, one must also study the sea surface topography in their 
immediate vicinity. 

As indicated above, the effect of age contrast between the lithosphere underlying young 
marginal basins and that underlying deeper ocean floor should also be included in a model of 
the regional geoid, Ideally, the geoidal undulations due to this effect should be superposed on the 
computed effect of slabs. However, sufficient data is not available to compute this age contrast 
effect for the entire 40^ X 40* region under study. In the preceding section, four marginal basins 
are modelled (see Figure 7). A number of basins in this region are unmodellcd; Figure 7 contains 
information gaps. It is therefore more useful to compute the geoidal contribution of these two 
effects (age contrast and anomalously dense slabs) in profiles along selected subsatellite tracks and 
then superpose them. This has been done. Comparisons with along-track sea surface profiles are 
shown in Figures 12a-12e. 

The ground track map of these five satellite passes is included in Figure 9. Figures 1 2a- 1 2c are 
Scasat 1 passes; Figures 12d and 12e areGeos-3 passes. The model geoids are quite smooth; they 
were computed at every third or fourth observation point. Observed data points are spaced, after 
smootliing, at 1 second intervals. Sea heights for the two Geos-3 passes are referred to an ellipsoid 
of semimajor axis 6378140 m and flattening 1/298.255. The Seasat passes are referred to an ellip- 
soid of semimajor axis 6378140 m and flattening 1 /298.2S7 An ocean tide model correction has 
been applied to Seasat heights but i.-ot to Gcos-3 heights. In computing model geoids in Figures 1 2a- 
12e, the same slab density anomalies (1.0 and 3.2 X 10® gm cm'^) used in Figure 9 were employed. 
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Figure 12c. See Figure caption 12a. Seahei^t profile is from the Seasat altimeter. 




Figure 12d. See Figure caption 12a. Sea height profile is from the Geos altimeter. 







Superposed on this slab effect was the effect due to age variation in the thermal lithospiierc (mar- 
ginal basin effect); this was computed using the very same model that was detailed in the preceding 
section (see Figure 7). 

The agreement in Figures 12a-l2e between computed and observed geoids is somewhat rough. 
However, in light of the major simplifying assumptions, a better agreement can hardly be expected, 
To obtain a good agreement the observed height profiles have been shifted downward by amounts 
from 4 to 10 meters. The need for this shift could be ascribed to a number Of things including: 

(i) compensation for orbit or altimeter bias error; (ii) unmodelled gravitational effects of very long 
wavelength; and (iii) compensation for a rather arbitrary choice of reference ellipsoid. A low degree 
and order background geoid could have been subtracted from the observed geoid, However, as dis- 
cussed below and above, this is not appropriate insofar as slabs appear to be capable of significant 
contributions to the very low degree-and-order field. If a hydrostatic (lower flattening) reference 
ellipsoid had been used, the shift needed to rectify theory and observation would be increased by 
about 30 to 40 m. Certain misfits can be readily interpreted as unmodelled effects centered in 
adjacent areas, For example, there is a long-wavelength low (. — 50 m) in the global geoid which is 
centered to the south of the region of study at a longitude of roughly 180° and a latitude of -70° 
This low could conceivably be associated with a component of tilt in the regional geoid of the 
southwest Pacific, In fact, two sea surface profiles, those of Figures 12b and 12e, have been ad- 
justed using a northwestward tilt of about 0,03 m/sec (~0,004 m/km). These two profiles are both 
ascending passes (see Figure 9). The other ascending pass, rev, no, 468 in Figure 12a, has not been 
adjusted and shows a sliglit tilt misfit with respect to the model result. 
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As indicated above, there are areas (marginal basins) which lack sufficient age data to conduct 
modelling of the effects of the local thermal lithosphere. This explains some of the misfit in the 
rip’'t halves of Figures 12c and 1 2d. 

A summary of slab density anomalies which provide an acceptable fit to observations is given 
in TaUe 2. The range (2.6 to 3.2 X 10® gm cm*^) in the New Hebrides slab density determination, 
is due to poorly constrained choice of density for the adjacent NG-NB-S slab complex. As indicated 
earlier, the 2.6 X 10® gm cm'^ density is slightly preferred. 

DISCUSSION 

This study has focused on interpreting a regional marine geoid of the southwest Pacific. There 
are some clear advantages to applying this same approach to a more global geoid. Estimates of slab 
density anomalies in this region, or any other individual region, should be somewhat improved by 
studying the entire marine geoid. Furtheniiore, a reasonable estimate of the global slab complex’s 
contribution to the geopotential power spectrum could be obtained. Recall that it was shown 
previously that slabs may, quite conceivably, contribute substantial power for degrees as low as 4 
and as high as 20 or greater. 


Table 2 

Slab Density Determination 


SLAB 

NOMINAL DENSITY ANOMALY 

Tonga-Kermadec 
New Hebrides 

1.0 X 10® gm cm’ 2 

*2.6 X 10® gm cm’^ 
**3.2 X 10® gm cm‘^ 


♦Adjacent NG-NB-S slabs lumped with New Hebrides in assigning densities. 
♦♦NG-NB-S slabs lumped witli Tonga-Kermadec. 
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As a test of this notion, a geoid was computed from a subset of GEM-9 gravity model coef- 
ficients (Lerch ct al, 1979) which consisted only of coefficients of degree and order 5 and below. 
This smoothed geoid retains most of the prominent 73 m high which is centered over New Guinea 
in the GEM-9 geoid map, The southeastern portion of this high was figured prominently in this 
study. In fact, this study indicates that a substantial portion of this high is due to dense slabs ‘ 
coupled with beind-arc elevation of the geothenn. 

It was also indicated previously that the New Hebrides slab density anomaly (Table 2) may not 
be typical of slabs worldwide. The nominal density anomaly of the Tonga-Kennadec slab is sub- 
stantially less (see Table 2). The South-Sandwich slab appears to have a density anomaly of no 
greater than 1.5 X 10® gm cm'® ; this estimate is obtained from a slab model fit to Benioff zone 
determination of Forsyth (1975) and Brett (1977). If the New Hebrides slab’s (and perhaps the 
NG-NB-S slabs’) nominal density anomaly is indeed quite higit with respect to most other slabs, 
then the slabs are not very substantial contributors to the global geopotential spectrum (degree 
< 20), This remains to be resolved, In any case, there obviously are important large scale gravity 
anomalies which are not associated with subduction zone density anomalies (see e.g., Cochran and 
Talwani, 1977, 1978 or McKenzie, 1977). 

Adopting the uniform density estitnates of 3.2 X 10® gm cm'® one finds that the undetached 
portion of the New Hebrides slab possesses an inferred negative buoyancy per unit length along 
strike of 1.1 X 10'® dyne cm''. Using 1.0 X 10® gm cm’®, one obtains a corresponding buoyancy of 
1 ,0 X 10*® dyne cm'* for the Tonga-Kermadec slab. Using the thermal model due to Molnar and 
Gray (1979) and McKenzie (1969), theoretical buoyancies of 1.2 X 10‘® and 2.8 X 10'® dyne cm’* 
are obtained for the New Hebrides and Tonga-Kermadec slabs respectively. It appears that the 
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Tonga-Kerniadec slab may be regionally compensated (cf. Griggs, 1972) to a large extent whereas 
the New Hebrides slab may not be. Recall that the New Hebrides slab has an exceptionally high 
sinking rate (convergence rate ~ 1 1 cm/yr according to Falvey, 1975: dip exceeds 60”). In esti- 
mating buoyancies from the thermal model, plate thicknesses of 100 km (Tonga) and 80 km (New 
Hebrides) were employed ; thermal and physical parameters were taken from Molnar and Gray ( 1 979), 

CONCLUSIONS 

Undulations in the altimetrically-obscrved geoid of the southwest Pacific are strongly controlled 
by (i) positive density anomalies in the subducting slabs of the region (see Table 2) and (ii) the 
effects of elevation of the geothemi in behind-arc lithosphere (corresponding to young marginal 
basins). Finer details of slab geometry do not obviously manifest themselves in the observed geoid; 
such gravitational effects are quite attenuated at sea level and are apparently mixed with crustal 
effects, oceanographic noise, etc. It appears that slabs in global composite may contribute sub- 
stantially to intermediate and long wavelength portions (down to spherical harmonic degree 3 or 4) 
of the earth’s gravity field. 
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